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Continuum and line emission of flares on red dwarf stars 

E. Morchenko 1 • K. Bychkov 1 • M. Livshits 2 


Abstract The emission spectrum has been calculated 
of a homogeneous pure hydrogen layer, which parame¬ 
ters are typical for a flare on a red dwarf. The ionization 
and excitation states were determined by the solution of 
steady-state equations taking into account the contin¬ 
uum and all discrete hydrogen levels. We consider the 
following elementary processes: electron-impact transi¬ 
tions, spontaneous and induced radiative transitions, 
and ionization by the bremsstrahlung and recombi¬ 
nation radiation of the layer itself. The Biberman- 
Holstein approximation was used to calculate the scat¬ 
tering of line radiation. Asymptotic formulae for the 
escape probability are obtained for a symmetric line 
profile taking into account the Stark and Doppler ef¬ 
fects. The approximation for the core of the H —a line 
by a gaussian curve has been substantiated. 

The spectral intensity of the continuous spectrum, 
the intensity of the lines of the Balmer series and the 
magnitude of the Balmer jump have been calculated. 
The conditions have been determined for which the 
Balmer jump and the emission line intensities above the 
continuum decrease to such low values that the emission 
spectrum can be assumed to be continuum as well as 
the conditions at which the emission spectrum becomes 
close to the blackbody. 

Keywords Red dwarfs, flares, Stark effect, H-a line 
core, escape probability, emission spectrum 
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1 Introduction 


Gcrshber g fe Pikcl’ncr) (1974) were the first who formu¬ 
lated the idea about the analogy of flares on red dwarf 
stars and on the Sun. At present it seems obvious that 
the physical nature o f the activity of th e Sun and red 
d warf stars is sim ilar ( Gershberel ( 2005 )1. 

Kostvuk fe Pikel’ner ( 1975 1 put forward the pioneer¬ 


ing idea about the important role of gas-dynamic effects 
in the formation of the emission of solar flares. The au¬ 
thors showed that the radiative cooling of the shocked 
gas can be a source of the flare’s emission in the op¬ 
tical range. The shock wave is produced by impulsive 
heating of the chromosphere by suprathermal particles. 

The Kostyuk-Pikel’ner model predicts the formation 
of a relatively homogeneous gas layer, a few kilometers 
of thickness, which is responsible for the flare’s optical 
emission in the continuum and spectral lines. The layer 
is located between the front of the downward shock and 
the region of plasma heated from above. 

This model was app lied to the problem of flares 
on red dwarf stars in IKatsova et alJ (19811) . The au¬ 
thors fulfilled calculations of a one-dimensional radia¬ 
tive shock wave and demonstrated the formation of a 
region with a small gradient of the plasma parame¬ 
ters: tem perature, pressure, an d density (see Figures 3 
and 4 in Katsova et al. (1981:)). This result is con¬ 
firmed by the calculations of the radiative cooling be- 
hind a shock fron t in th e at mospheres of c o ol sta rs (see 
Fadeyev fe Gillet ( 2001 ) and Belova et al. ( 2014 )) with 
a more detailed account for the elementary processes in 
the plasma: ionization, recombination, bremsstrahlung, 
excitation and de-excitation of discrete levels of atoms 
and ions as well as radiation scattering at the frequen¬ 
cies of spectral lines. 

Red dwarf stars have a complex structure of flare 
spectra: the intensities of the emission lines strongly 
increase, the continuum emission is also enhanced; it 
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is well approximated with Planck curves obtained for 
different temperatures at the blue and red si d es rel a- 
tive to the Balmer jump (see iKowalski et al. C 2013 h). 
Layers producing the blackbody emission should be in 
a state close to thermodynamic equilibrium, whereas 
the gas generating strong emission lines is essentially 
non-equilibrium. Thus, the flare plasma includes at the 
same time equilibrium and non-equilibrium regions. 

Therefore, we investigate the problem of the forma¬ 
tion of an emission spectrum in a wide range of the 
cloud parameters. We consider a homogeneous gas 
layer of thickness C with definite values of temperature 
and density and consisting o f pure hydrogen. 

We develop the model of Drake fc Ulrich ( 198Clh in 
three points. Firstly, we consider the influence of the 
layer’s radiation (bremsstrahlung and recombination) 
on the occupation of atomic levels. This development 
is necessary, as the flare luminosity is stronger in the 
optical range than that of the quiescent atmosphere of 
the whole star. Secondly, we take into account that 
broadening of the high members of spectral series is 
caused exclusively by the linear Stark effect. Thirdly, 
we consider the two-temperature approximation for a 
flare layer. 

In Section 1 the problem is formulated. In Sec¬ 
tion 2 we write the set of equations of ionization bal¬ 
ance with regard to the continuum. In Section 3 the 
calculation of the continuum emission is described. In 
Section 4 we present the method of calculation of the 
quantum escape probability. In Section 5 the intensities 
of the emission spectral lines are calculated. Section 6 
presents the results of the emission spectrum calcula¬ 
tions of a layer with parameters possible in the flare 
regions of red dwarf stars. 


2 Statement of the Problem 

We calculate the emission in spectral lines and contin¬ 
uum of an homogeneous pure hydrogen layer. Its tem¬ 
perature, density and dimensions are assumed to be 
given. The value of the ion-atom temperature, T a i , is 
taken to be greater than the value of the electron tem¬ 
perature, T e , because both ions and atoms are heated 
mor e intens i vely tha n elec t rons at the shock front (see, 


e.g.. IZcrdovich fc Raizcr (119671 11. 


We use the results of the radiative shock wave cal¬ 


culations bv lFadevev &; Gilletl (1200111 and I Belova et al 
(120141 to state the input parameter set. These calcu¬ 
lations were fulfilled under the conditions of the atmo¬ 
sphere of a red giant star with an effective temperature 
of To = 3000 K. The values of T a i and T e at the shock 
front are respectively 7 — 10 eV and « 0.7 eV. 


The Lyman-a (Ly-a) optic al depth, Tys , in the cool¬ 
ing region, as it follows from Belova et alJ (20141. lies 
in the range 


1T0 7 < ri2 < 3-10 7 . 


(1) 


The value of T 12 weakly depends on the density of the 
gas, as the cooling rate, dS/dt , is proportional to the 
hydrogen atom concentration, N a . Thus, N a is can¬ 
celled: 

dt dS 

g?ti 2 oc N a u— dS, —— oc N a => dr \2 oc udS . (2) 

dS dt 

Here, 


S = T a . 


Ne_ 

N n 


-T, 


( 3 ) 


N e is the electron density, u is the velocity of the gas 
flowing from the shock front, and dt is some time inter¬ 
val from the moment when a given gas layer intersects 
the shock front. 

The value of TVh is equal to the sumrna of the proton 
and atom concentrations: 

N u = N p + N a . (4) 

We accept Nr typical for the M-dwarf chromospheres: 


3-10 14 cm ** ^ Nn ^ 3-10 1D cm 


.-3 


16 


-3 


( 5 ) 


The layer thickness, C, is defined to be such that the 
condition ID is satisfied. The electron density is equal 
to 


N e = N p . 


( 6 ) 


The values of T a i and T e were obtained from a red 
giant star, which has a lOOx smaller gravity compared 
to a red dwarf star. The ion-atom temperature at the 
shock front is determined by the shock velocity, uq, and 
weakly depends on the density of the gas. If the value 
of uq lies in the range 50—70 km/sec, then it is true 
that 


6 eV ^ T n ,i ^12 eV. 


( 7 ) 


The value of T e at the shock front is determined by 
the temperature Tq of the unperturbed gas and heating 
due to the compression. The evolution of T e is caused 
by the balance between the supply of energy in elas¬ 
tic collisions with ions and atoms and radiative losses 
in the Lyman-a line. The radiative cooling in the hy- 
d rogen lines is the most intensive at T e ~ 1 eV (see 
Belova et al. ( 2014 )). Therefore, we assume that 


0.8 eV <T e < 1.5 eV. 


( 8 ) 
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3 Balance Equations for the Populations of the 
Levels 


The following processes we take into account: the elec¬ 
tron impact ionization, excitation, and de-excitation, 
the triple recombination, radiative recombination (in¬ 
duced and spontaneous) as well as the influence of the 
layer’s bremsstrahlung and recombination radiation on 
its ionization and excitation states. 

Th e impact coefficients are taken from Ijohnson 
( 1972j). We use t he Kramers approximation (see, e.g., 
Sobel’man ( 19631) 1 for the photoionization cross sec¬ 
tions from the excited levels. The rates of the inverse 
processes are calculated following the detailed balance 
principle. The spectral lines optical depth can be rather 
large in this problem, so the scattering of line radiation 
should be taken into account. We introduce so called 
effective Einstein transition coefficients for the transi¬ 
tions between hydrogen levels i and k by the following 
formulae: 


f A ki 1 

| 

r A k i i 

B *ki 

= 

B ki 

{ B* k J 

1 1 

{ Bi k ) 


x Oik ■ 


( 9 ) 


Here, 0i k is a photon escape probability out of the pro¬ 
cess of scattering. The values of () lk are calculated in 
the section [5] for the resonance as well as subordinate 
transitions. 

According to the standard notation, A k i and B k i 
are respectively the probabilities of spontaneous and 
induced transition from the upper level k to the lower 
level i, and Bi k is the absorption coefficient. 

We write the set of the balance equations for the 
discrete levels of a hydrogen atom: 

k— 1 n 

Y, aijNj + a lk N k + 53 aijNj = aiN e N p , (10a) 

3 =1 j=k +1 

k— 1 n 

53 CL 2 jNj + a 2k N k + 53 a 2 jNj = a 2 N e N p , (10b) 

3 =1 j=k+1 


k— 1 n 

^ ^ Q'kjNj H- a kk N k ^ ^ a kj Nj — akN e N p , (10c) 
j =i j=k+1 


k— 1 n 

y a n j Nj + a nk N k + y 

&nj N 0 = a n N e N p , (lOd) 

3 =1 j=k+l 


We assume that the maximum value of the principal 
quantum number, n, is equal to an effective series li mit. 
It can be evaluated by the Inglis- Teller formula Vidal 
( 196filW as 


lgn = 3.10 — 0.13lg(2V e /cm 3 ). 


( 12 ) 


The following designations are introduced: 

N k — the density number of the hydrogen atoms on 
the discrete level with the principal quantum number 

k; 

r k , 7 fc — the radiative and triple recombination coeffi¬ 
cients on the fc-th level, respectively; 
dik (i ^ fc) — the transition rate from the level k to 
the level i; 

Q>kk — the escape rate from the level k. 

We write the explicit formulae for the matrix coeffi¬ 
cients in the left side of the equation set. 

Excitation rate ( k<i ) is: 

a ik = ~ [q ki (T e )N e + 4nI ki BlJ[ , (13) 


where 

Qki (T e ) is the electron impact excitation; 

I k i is the spectral intensity of the layer continuous ra¬ 
diation in the frequency of fc-f-H transition. 
Deactivation rate (k>i) is: 


a ik = - [q ki (T e )N e + A ki + 4tt I ki B * ki \, (14) 

where q k i(T e ) is the impact deactivation coefficient. 
The diagonal elements are 

k- 1 

a kk - ^ [q kj (T e )N e + A* kj + 4t Tl kj B* kj \ + 

7=1 

n 

+ 53 [q kj (: T e )N e + 4t Tl kj B* kj ] + V k + q k ■ N e , (15) 

j=k+1 

where q k is the electron impact ionization coefficient, 
and V k is the photoionization rate from the level k: 

r°° t 

V k = 47t / y— <y k {y) du . (16) 

Ju k 

Here, v k is the photoelectric threshold from level k, h 
is the Plank constant, <j k {v) is the cross section in the 
Kramers approximation: 

1 <17) 


where 

oti = n + 7 ifV e , i = l,...,n. 


where a is the fine structure constant, ao is the Bohr 
radius. 

(11) The spontaneous recombination coefficient, , can be 
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calculated from tm> using the principle of detailed bal¬ 
ance: 


rf =npk)Er{Pk), /E = Jfr; 

(18a) 

lZ{x) = C r ■ x 3 / 2 e x , C r = ^— </E a 4 cag , 

O y O 

(18b) 

where E\{x ) is the exponential integral of order 1. 

The coefficients of induced recombination, r[ md ' 1 , is 
given in the next section. It is true that 

(s) (ind) 

rk = n' + r K k 

(19) 

Summing the set (flOl) over all discrete levels yields the 
identity 

E + * N *) • N i = E ( r ‘ + 7iJVe) • 

(20) 

As a result, we come to the set of homogeneous equa¬ 
tions with singular matrix for the relative level occupa¬ 
tions, yi=Ni/N a : 

n 

E ^ijVj = °> * = !,•••,«• 

1=1 

(21) 

where 


&ik — (Lik (J^k H - Qk-^e) 5 

(22) 

n 

ai = n+ liNe, T = E “!■ 

1=1 

(23) 

To avoid the singularity, the last equation (fc= 
placed by the normalization condition 

-n) is re- 

T-1 

II 

-Sfi 

(24) 


fc =1 


We take into account the scattering of the ionizing ra¬ 
diation using the division of r k by 1 + Tfe, where r k is 
the optical depth on the photoionization threshold from 
the level k. 


4 Continuum Emission 

Continuum emission consists of two components: free- 
free and bound-free. The following formula is valid for 
an homogeneous layer: 

E = S v • (1 - e- T ") , (25) 

where S u is the source function, which is equal to the 
ratio of the emission coefficient, £„, to the absorption 


coefficient, x v \ 


S v = — = 


+ EE 

j = k 


lb) 


E • 

j=k 


„(0 


(26) 


The summation is over all the levels for which the fre¬ 
quency of the photoelectric threshold is less than v: 


Ry 

The Menzel factor is designated by b k '- 
N k (2tt m e T e f /2 


bk = 


N e N p 


k 2 h 3 


exp (-/?*). 


(27) 


(28) 


It shows the deviation of the level population from its 
equilibrium value. This factor is equal to 1 when exci¬ 
tation and ionization are described by the formulae of 
Boltzmann and Saha. 

We use the K ramers app r oxima tion for the absorp¬ 
tion coefficients ( Sobel’man (1963), §34.4): 


= 


2 8 

3-\/3 


(i)_ 2 9 5/2 51, 

'Afe 3y^3 a7F a 0 ( ’ I 1 m ) 




hv 


3/2 


1 


1 - — e~ y ) e fih N e N t 


bk 




T e J 

y > Pk, 


(29a) 


(29b) 


hv 

where y = — . 

-L e 

We write the relation between the ionization coef¬ 
ficient from the level fc, , and the bremsstrahlung 
coefficient, k v 


lb). 

IS • 

b k - e~ v 
Ik) ~ ^ k i _ e -v ’ 


>k 


<Pk = -^P kePk - 


(30) 


The value of t v is equal to the sum of bremsstrahlung 
optical depth, rl, and photoionization one, tE : 


t -r (&) 

‘IS 1 u 


t E = I K 


■(b) 


E- 

j=k 


ii) 


■C. 


(31) 


The free-free absorption and emission coefficients are 
related by Kirchhoff low: 


£ (6) 

-W = B ^- 

ytv 


( 32 ) 
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The Menzel factor connects the radiative recombination 

(r) 

coefficient, £ vk , and photoionization coefficients for the 
fc-th level: 


(r 


(i) 


B v {T e ) 


l~e~ y 
b k - e~y ’ 


(■y > Pk) ■ 


(33) 


Here, u is the angular frequency, a u is the spectral line 
profile, and k u is the line absorption coefficient. R 
denotes the radius-vector of a point belonging to the 
surface S. 

Let the gas occupies a layer of finite thickness C. In 
this case, the mean escape probability is equal to 


Using the last formula and taking into account (EZl) , we 
can transform (|26D : 


6 = j du a u G(T u ), 


(37) 


S v = 


£ ] /xi b) + £ /*£,•) • 

l+Ex«/x« 

j=k 


= B u (T e ) 


jhc xl b) b j ~ e~ y 

n >£) 

1 +y^~ 

A, (b) 

J = k 


= By(T e ) ■ ip k (y), 


(34) 


where the following function ip k is introduced: 


V'fc(y) = 


1 + D Tj 

3—k 


1 + 1 _ e -y s <Pj ■ ( b j - e y) 


(35) 


j=k 


The induced recombination coefficient is equal to 

-T’ = K(ft) ■ r £<%• (36) 


/3fc 


— 1 y 


As is easily seen from (0H1) . the continuum optical 
depth, r„, is the function of y when the frequency is 
varied and the other parameters are unchanged. 


5 Photon escape probability 

The source of optical emission moves slowly and co¬ 
hesively, i.e. there is no noticeable velocity gradi¬ 
ent Ther efore, the Biberman-Holstein approximation 
( Holstein ( 1947h l is used to calculate the scattering of 
line radiation. Consider a stationary gas that occupies 
a volume V bounded by a closed surface S. Let r is the 
ra dius-vector of a p o int wi thin the volume. According 
to iBibcrman ct al. (|1982lf , the photon escape proba¬ 
bility from this point outside the homogeneous plasma 
can be written in the form: 


where 


Q{r u ) = —[l-2Ek{2r u )\. 


ku,C 


Tcj = 


(38) 


In these expressions, t w is the optical depth at the cen¬ 
ter of the layer an d E%(x) is the th ird-order exponential 

integral (see, e g.. ISmirnovI (11972 )). _ 

As stated in Sobolev &; Grinin ( 1996h . the very wide 
wings of the Balmer lines in the optically thick flare 
plasma are due to the linear Stark effect. When the 
Stark broadening by electrons and ions is treated in 
a semi-classical picture, the static li mit is vali d i n th e 
distant long-wavelength line wing ( Sobel’man ( 19631) . 
§38.1): 


afpdu ■ 


15 (B kn ,£ 0 y/ 2 
4\/27r Aw' 


5/2 

kn' 


du>. 


(39) 


Here, B kn ' is the Stark broadening parameter of the 
transition from level k to level n' (k > n'), £q is the 
normal Holtsmark field, and Aui kn i = u> — ui kn i, u> k ri 
is the transition angular frequency. The Stark width is 
approximately equal to 




(40) 


The spectral line profile is determined by radiation 
damping, Doppler effect and pressure effects. The ex¬ 
tensive tables of the Stark broadened profiles for hydro¬ 
gen pl asma at T a , = T e have been published bv IStehlc 
(1994:). In the case of T a i > T e we will use a simple 
model of the line profile with the Doppler core and the 
Holtsmark wings. 

We consider, e.g., the H-a line profile in order to 
justify this model. The following set of parameters is 
assumed: 


14 


T a i = 6 eV, T e = leV. 


(41) 


/• f (dS, R r) 

<9(r) = J dio j) 47r |R„ r |3 exp - r D • 


We introduce the value of f2 32 such that 

-l 


H 32 = v 


e 2 a 0 


, where v is the hydrogen atom 
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velocity (along a line of sight). The value of v can be 
evaluated as 


'3 T ai 
m H 


n 32 ~ 2.87- 10 12 sec" 1 . 


(42) 


Here, mn is the hydrogen atom mass. The Doppler 
width, A u>k n ,, is equal to 


Aw fen / OJ k r 


Here, k = 3, 
,d 


2 T a . 


itihc- 


.2 ‘ 


(43) 


= 2 => Au£ 


32 


3.246 • 10 11 sec“ 


and A oj ^ 2 ^ ^ 32 - Therefore, the static limit (1501) can 
be used wh en AW 32 ( 2 > ^ 32 , i.e., in the line wing (see 
Sobel’manl Jl963l) . §38.1, §36.6). 

On the other hand, we introduce the parameter 7 , 
which is equal to the sum of the radiative damping con¬ 
stant, 7 rad, and the parameter 72 : 

7 = Trad + 72 ; 72 ~ 72 tt 3 N V - 1 ~ 5.4 • 10“ 4 7V, (44) 

where N is the concentration of perturbing particles. 
Trad < 2 • 4.67 • 10 8 sec " 1 => 7 ~ 72 . If N ~ N e then 
7 /(A-C 1. Therefore, the central part of the I l-o 
line profile is the Doppler profile (see lSobcl’man ( 19631) . 
§36.6). 

In general, there is some analogy between the line 
profile with the Doppler core and the Holtsmark wings 
and the Voigt profile when the Voigt parameter a is 
much less than 1. The normalized Voigt function is 


aY = 


1 


V^Awg,, 


- H(a , v) 


(45) 


where v is the dimensionless frequency displacement 
from the line center in units of Aw® n ,. For a <C 1 it 
can be assumed t hat in the first approximation (see, 
e.g., Ivanov ( 19691 ) I 


J exp (—f 2 ), 0<M<|vo| 

H(a,v) = l a n | | 

\-j= M>M- 

t \rkv z 

The value of |uq| is the solution of the equation 


(46) 


5.1 Ly-a line 

The following model symmetric profile is similar to the 
Voigt profile (j45]l : 


C 


v^Aw®, 


®{bkn',v) , 


(48) 


where C is a normalizing constant. The function $ is 
defined by 

f exp (—v 2 ), 0 <H<|u o | 

$(bkn',v) = < 15 b kn , 

[ 1 T 2 W - N - 

In these expressions, the value of parameter bkn 1 is sig¬ 
nificantly less than unity: 

bkn^(^) «1. (49) 


We can find the constant C from the equations 


+ OO 


1 = / a u dui = C ■ erf vq + 


5 b 


'kn' 


vl /2 ) ’ 


15 bkn' 


4V2 v 5 / 2 


= exp(-vfi). 


(50) 


(51) 


iDrake & Ulrich (1980) introduced the model profile 
phenomenologically. Here, we give the theoretical foun¬ 
dation for its application. 

To make an estimate of the optical depth in the Ly-a 
line at |u| = |i>o|, we use (UTl) . According to (HOI) . the 
Stark width of the Ly-a line is approximately equal to 
& 2 i£o = 01-10 10 sec" 1 , the Doppler width is Aw^ = 

1.75 • 10 12 sec" 1 ; whence it follows that 

b 21 = 4.385 • 10 " 4 « 1. 

Thus, the approximation (H51) is adequate in our prob¬ 
lem. Using (1501) and (1511) . we find that |uo| ~ 3.095 and 
C = 1. If the optical depth at the line center is equal 
to t/O = 3 • 10 7 , then it is true that 


exp (-Uq) = 


V^o 


(47) 


The transition from the Doppler core to the Lorentz 
wings in the Voigt profile is quite sharp (iMihalas 
( 19741) 1 . 


'Tw(M) = t U0 • exp(—u 0 ) = 2 • 10 d » 1 . 

Therefore, the Ly-a photon lives the medium at the 
frequencies of the distant line wings. 

As a result, equation (1371) can be written in the form: 


+ OO 

15 f dv 

e i2 =b 2 l-/= J G{Tu)^j 2 -, T UJ =T ul {v). 


( 52 ) 
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It is obvious that 


Then the principal quantum numbers, k and n ', must 
satisfy the condition 


, d 15 & 21 

4^2 ' v 5 / 2 ’ 


n 


15 


^21 


4\/2 w 5 / 2 


Here, is the absorption coefficient at the Doppler 
profile center: 


n 2n 2 e 2 

=- fn'kN, 


1 


m e c ' ^/irAuj]? n , ' u ° 2 


n 

= ( 53 ) 


/„/fc is the oscillator strength for absorption (n' = 1 , 
k = 2). Taking into account that 1, we obtain 

the asymptotic formula for the mean escape probability: 


f.2/5 


\4V5/ V^(r£) 3/S 




2/5 ' 

tJ 


Consider the integral 

+oo +oo 

/' dr 1 r dr 

‘= / ew- =i y [l-2B 3 (2r)|—, (0,1). 

0 0 

By using the formula of integration by parts, this equa¬ 
tion can be represented in the form: 


+oo 

If dr 

‘--J Wr)- = - 


dr 1 [1 - 2£ 3 (2t)] 


4/t 


+oo 


If the value of k lies between 0 and 1, then the r.h.s. of 
the last equation tends to zero. Thus, we obtain that 


r /x 2 r(i — i/ n) ______ 

2/i-l 2 1-1 /^ ’ ^~ K ' 


(54) 


where T is the gamma function. In this problem, the 
value of n is equal to 5/2. Finally, we come to 


3/5 


1 


712 


( b 21 £ (A _ 

V Aw 2 D ! J (T°)W 


(55) 


The exponent 3/5 in the last equation is the feature of 
the linear Stark effect. 

5.2 Generalization for the case of a multilevel atom 

Strictly speaking, the resulting expression for 0x 2 is 
valid only for a model two-level atom. However, the 
escape probability approach can be applied also for a 
multilevel atom. We will use the mode l profile (l48l) up 
to b = 3 • 10 -2 (for more details, see Stehle Jl994l ff. 


3/2 


b > bkn' = 


&21 ■ 


(56) 


Constant C is approximately unity in the range 10 4 ^ 
bkn 1 ^ 10 -2 ; therefore, it is true that 


3/5 


1 


/ Bkn '^0 

{ A^ri 


(57) 


where r// is the optical depth in the n' k transition. 

Since the Stark width increases with increasing k, the 
Stark broadening dominates over the Doppler broaden¬ 
ing in the limit k 1. Then we use the following model 
profile: 



0 < 


< 


2 05/2 ’ 


>|A)| 


(58) 


UJ — 0/0 

Here, ao is some constant, (3 = —-—. The normal- 

Bkn'C'O 

ization condition gives 

0.1 

A) ~ ^ 5 ~ "75 ■ 

Bkn'C' 0 

After some calculations we obtain that 
0.36 


7 r 2 e 2 


^ n'k ~ /r 5 ^cjn — fn'kNn'&uo^' • 

r d / 5 m e c 

7 ^o 


(59) 


5.3 Application to the observations of the H-a line 
profile 

It is interesting to note that the profile with the Doppler 
core and the Stark-broadened wings is similar to the H- 
a line profile in the flar e on UV Ceti (dM5.6e) published 


Eason et al . (1199211 . The authors found a Doppler 


profile to fit the line core: for H-a during outburst, the 
Doppler width was equal to 0.9 A. Thus, it is true that 
T ai w 8.8 eV and A uj 32 « 3.93 ■ 10 11 sec -1 . Further, 
Eason et al. found that log (N e ) = 14.75 => B 32 £ o ~ 
5.33 • 10 10 sec -1 and 7 ~ 2.5 • 10 11 sec -1 . Thus, we 
obtain that 7 /Aa;^ ~ 0.6 and 632 -C 1 . 

The authors satisfactorily represented part of the 
long-w avele ngth wing w ith a Stark profile from the pa¬ 
per by 1 Underhill k , Waddell (11959 1. In this paper the 
Stark profile was calculated taking into account only 
the ion broadening in context of the Holtsmark dis¬ 
tribution of the electric microfields, while the Stark 
broadening by electrons was neglected. But as stated in 






































Sobolev fe Grininl (119991 . the electron broadening dom¬ 


inates in stellar flares. We consider that the deviation 
of the observed Ha line profil e from the St a rk pr ofile 
in the red wing (see Fig. 9 in Eason ct al. (|1992 1) is 
almost entirely due to neglecting the contribution of 
electrons to the line broadening in the optically thick 
flare plasma. 


6 Spectral Line Emission 


To simplify the notation, we have in mind the Balmer 
series. The intensity of the emission in the spectral line 
frequencies, ll k \ is equal to 



(60) 


where index k identifies the upper level of the transi- 
tion k—>2. We have denoted by r,) the optical depth 
in the frequencies of the 2 —>k transition; the negative 
absorption being taken into account: 

T^ fe) = ■ a v - C. (61) 


Here, 


(fc) _ Tre 

— -/ 2 fc ‘ 

m e c 



(62) 


The spectral line intensity, I k , is equal to the integral 
of (ITOl) over the line profile: 


line emission merges into continuum. They are defined 
by the following condition: 

=T < f ) {v u ) =T v (v k ). (66) 

The simplified formulae (l65ll and (IMll do not introduce 
significant inaccuracies in Ik , as the transition region is 
quite narrow. Moreover, the emission intensity in this 
region is much lower than the intensity in the line core. 

Note that there are some factors simplifying the cal¬ 
culation of integral (lC)5l) . If the condition (P) is satis¬ 
fied, then the central optical depth, Tfc. is rather large. 
Consequently, the intensity in the line core frequencies 
coincides with a source function. 

In those parts of the line profile where is not too 
large, we use for a„ a simple approximation (1581) . as¬ 
suming that for the frequency displacement up to /?o, 
the layer remains opaque. We introduce frequency dis¬ 
placement /3 0 > 3 which satisfies the condition 

t-W > 5. (67) 

In the offset range 

0<|/8|<|A)|, (68) 

it is assumed that the intensity is equal to the source 
function, and in the line wings we calculate the integral 

r@u 

D k = 2/ (l - e~ T ^ Uk) ■ e~ TkrW2 ) d^. (69) 

Jpo V ' 


Ik 



1 



dv. 


(63) 


In this section, we consider that the continuum absorp¬ 
tion is negligible within the profile of a spectral line. 
So, we write the expression for the source function tak¬ 
ing into account only the line emission and absorption: 


Here, parameter r k is defined as 


(fc) 

Tfc = ^o 


3 27T 

2 ' B k2 £ 0 ' ’ 


(70) 


and the upper limit is calculated according to (1S51) . 
Thus, the total energy emitted in the line is 


SP=B v (T e )- ) 6 V (64) 

h/bk - e y 

The Planck function varies very slowly within the line 
profile. Therefore, the source function at the frequency 
of the line center, Sk, can be taken outside the integral 
sign: 


I k = S k -(2f3 0 +D k )-^—, (71) 

Z 7T 

where the Stark width, B k 2 £ o, is determined by formula 

(TO . 


7 Results 


h 



vi 


(65) 


For the same reason, we consider that the continuum 
optical depth within the line profile is equal to T v (yk). 
The limits of integration are the frequencies where the 


We performed calculations of the spectral line and con¬ 
tinuum intensity for different values of density and tem¬ 
perature of a homogeneous gas layer under the condi¬ 
tions of flares on red dwarf stars. The layer thickness C 
is fixed. We calculated the optical depth at the center 
of the Balmer series lines (see table [T]) and the Menzel 
factors; in all cases the values of b k do not differ from 
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unity (see table [2]). It is important that there is no ex¬ 
ternal source of radiation ionizing the hydrogen atom 
from the excited levels. 

Let us address the results of calculations of the emis¬ 
sion spectrum. Figure 1 shows how the emission lines 
are “immersed” in a continuum with increasing density 
of the gas. The horizontal axis represents the wave¬ 
length A in Angstroms, vertical — the logarithm of the 
ratio I v /B v . Balmer lines are indicated by vertical line 
segments, its height is proportional to the line center 
intensity. In all the graphs, the following values of the 
parameters are assumed: T e = 1.2 eV, T ai = 6 eV, and 
C = 20 km. The intensity of the emission at the centers 
of the Balmer series lines is very close to the Planck 
function at the electron temperature B\(T e ). There¬ 
fore, all the spectral lines have the same upper bound¬ 
ary corresponding to zero on the vertical axis. As the 
rise of the continuum occurs, first H a “immerses”, then 
H /3 and finally high members of the Balmer series. This 
is due to the fact that continuum absorption coefficient 
increases toward longer wavelengths. 

We temporarily exclude spectral lines in order to 
more conveniently trace how the continuous spectrum 
depends on the layer parameters. Figure 2 represents 
the ratio I\/B\ for two values of the electron temper¬ 
ature and three values of density. Atomic-ion temper¬ 
ature is equal here to 10 eV, and the thickness of the 
layer — to 10 km. The values of the optical depth in 
a continuous spectrum at the wavelength 4170 A are 
shown on each curve. It is seen that increasing of the 
continuum opacity essentially depends on the electron 
temperature: at higher values of T e , the transition to 
the Planck curve occurs at lower values of TVh. This fact 
is caused by increasing of the occupation of excited lev¬ 
els (fc > 3) with increasing the electron temperature. 
Indeed, the value of oc TV*,, as it follows from equa¬ 
tions (l28l) and (129bl) . 

While the gas density is increased, the intensity of 
the continuous spectrum I\ tends to the Planck func¬ 
tion, remaining less than it. The curve portions that 
are closer to the Planck curve are located in the red re¬ 
gion of the spectrum and at the blue side of the Balmer 
jump. When the density increases further, the blue 
part of the spectrum (redder at the Balmer jump) ap¬ 
proaches the Planck function. 

The shape of the line core is largely determined 
by the ion-atom temperature, as the Doppler width, 
Awj^,, depends only on the square root of T ai (if the 
principal quantum numbers fc and n! are unchanged). 
It is true that at higher values of T a i, the approxima¬ 
tion for the core by a gaussian curve is valid for a higher 
number of members of the Balmer series. 



Fig. 1 The continuous spectrum and the height of the 
lines of the Balmer series for the four density values. 
From bottom up: 1-10 15 cm -3 , 3-10 15 cm” 3 , 1-10 16 cm” 3 , 
3-10 16 cm” 3 . Under each curve its number is written. The 
upper curve almost coincides with the Planck function; the 
continuum absorption is not negligible in this case. 



Fig. 2 The continuous spectrum normalized to the Planck 
function (spectral lines are excluded). The dashed lines cor¬ 
respond to T e = 0.8 eV, solid ones — to T e = 1.0 eV, 
Numbers I, II, and III correspond to the following values 
of the total number densities of the hydrogen particles: 
2 • 10 15 cm” 3 , 7 ■ 10 15 cm” 3 , and 3 • 10 16 cm” 3 , respec¬ 
tively. Next, the values of the optical depth at wavelength 
A = 4170 A are shown for each curve. 

Table 1 Decimal logarithms of the optical depth at the 
centers of the Balmer lines (T e = 1.2 eV, T a i = 6 eV, and 
£ = 20 km) 


log A 

H a 


h 7 

H 5 

15.0 

6.40 

5.54 

5.06 

4.73 

15.5 

6.90 

6.04 

5.56 

5.23 

16.0 

7.40 

6.54 

6.07 

5.74 

16.5 

7.90 

7.04 

6.57 

6.24 
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8 Discussion 


Acknowledgments 


The main result of our work is the calculated plasma 
emission spectrum in a broad interval of the contin¬ 
uum optical depths, from the transparent gas to the 
gas whose emission intensity is close to the Planck func¬ 
tion. The proximity to the Planck function in the case 
of strong self-absorption is due to the mechanism of im¬ 
pact ionization and excitation (from electron collisions) 
of the plasma, whereas radiative processes are sec¬ 
ondary. The approach used in this paper fits well for red 
dwarf stars with their effective temperature of about 
3000 K. However, the case of G- and F-type dwarfs 
requires taking into account the photospheric emission, 
e.g., photoionization of hydrogen atoms from exited lev¬ 
els. In the case of appreciable emission and absorp¬ 
tion in the continuum its (continuum) influence must 
be taken into account in the calculation of spectral line 
int ensities. Here o u r wor k dev elops the approa c h use d 


Drake fe Ulrich ( 198Clh and Katsova et al. ( 1991 ). 


where the line intensities were calculated under the as¬ 
sumption of the transparent continuum bringing small 
contribution to the emission. 

In the framework of the considered problem the in¬ 
tensity of the subordinate lines is close to Planck for 
any layer parameters implemented during the flares on 
red dwarf stars. Therefore, the relative height of the 
lines depends on the optical depth in the continuum. 
If it is small then strong lines are observed against the 
background of faint continuum, which is described by 
formulae valid for transparent gas. In the case of appre¬ 
ciable self-absorption in the continuum its spectrum be¬ 
comes close to Planck; thereby, the lines weaken. How¬ 
ever, in the homogeneous layer model no strong lines 
on the background of th e Planck continuum observed 
by Kowalski et al. ( 2013h are obtained. Probably, we 
see in the flares the total emission from several layers; 
some of them are transparent in the continuum, while 
other ones are not. We hope to get a better agreement 
with observations by calculations of the radiative cool¬ 
ing behind the shock front under the conditions of the 
red dwarf chromospheres, since the shocked gas can in¬ 
clude at the same time equilibrium and non-equilibrium 
regions. 


E.S. Morchenko is grateful to Prof. Yu.K. Zemtsov 
for valuable remarks. This work was partly supported 
by the Russian Foundation for Basic Research (project 
code 14-02-00922), and the Scientific School (project 
code 1675.2014.2 NSh). We thank the anonymous ref¬ 
eree for very useful suggestions and comments. 


Table 2 Menzel factors for the discrete levels of the hy¬ 
drogen atom 


log A 

1 

2 

3 

4 

5 

15.0 

1.030 

1.031 

1.006 

1.001 

1.000 

15.5 

1.009 

1.009 

1.003 

1.002 

1.001 

16.0 

1.003 

1.003 

1.002 

1.001 

1.000 

16.5 

1.001 

1.001 

1.001 

1.001 

1.000 
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